Angular-resolved photon-coincidence measurements in a multiple-scattering medium 
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We present angular-resolved correlation measurements between photons after propagation through 
a three-dimensional disordered medium. The multiple scattering process induces photon correlations 
that are directly measured for light sources with different photon statistics. We find that multiple 
scattered photons between different angular directions with angles much larger than the average 
speckle width are strongly correlated. The time dependence of the angular photon correlation 
function is investigated and the coherence time of the light source is determined. Our results are 
found to be in excellent agreement with the continuous mode quantum theory of multiple scattering 
of light. The presented experimental technique is essential in order to study quantum phenomena 
in multiple scattering random media such as quantum interference and quantum entanglement of 
photons. 
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When coherent light propagates through a randomly 
disordered medium such as biological tissues, turbid me- 
dia or powders, the light waves are multiple scattered and 
interfere [1, 2]. The hallmark of multiple scattering is the 
formation of a complex intensity speckle pattern in the 
light transmitted and reflected through the medium. Al- 
though these speckles exhibit large intensity fluctuations 
and appear to be random in space and frequency, these 
fluctuations can be statistically correlated. The resultant 
intensity correlations provide important information on 
the wave transport parameters such as the diffusion con- 
stant of light [3] . In the diffusive regime of multiple scat- 
tering, the light transport can be described by a diffusion 
equation and interference effects vanish after ensemble 
averaging over all realizations of disorder. The ensem- 
ble average can be experimentally realized by studying 
many independent speckle patterns corresponding to dif- 
ferent sample configurations. However, in very strongly 
multiple-scattering media, deviations from diffusive light 
transport appear leading to non- vanishing mesoscopic in- 
tensity correlations and ultimately to Anderson localiza- 
tion. These effects can be explained by classical wave 
interference, neglecting the quantum nature of light. 

An intriguing question both for reasons of fundamen- 
tal interest and potential applications is how multiple 
scattering affects the quantum properties of light. Trig- 
gered by the theoretical quantum optical framework of 
multiple scattering [4], many proposals have been put 
forward addressing new phenomena that are of purely 
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quantum origin. In particular, the contributions from 
quantum vacuum fluctuations have to be considered in a 
quantum optical description. It has been proposed that 
spatial quantum correlations can be induced by multi- 
ple scattering [5] arising between two separated spatial 
or angular optical modes in the far field of the multiple 
scattered light. By extending this theory to several in- 
cident quantum states, it was predicted that quantum 
interference may be induced by multiple scattering [6-8] 
and even survive the process of ensemble averaging [9] , 
which is of potential interest for quantum information 
processing. Recent theoretical work demonstrated the 
use of multiple scattering to enhance the information ca- 
pacity for optical communication by exploring the intrin- 
sic coupling between input modes and output modes [10]. 
This effect can be even further increased exploiting the 
quantum optical aspects of multiple scattering using non- 
classical light sources that exhibit sub-Poissonian photon 
statistics [11, 12]. 

The quantum optical properties of multiple scattered 
photons can be accessed experimentally through their 
statistics and correlations [12-16]. So far, only ini- 
tial experimental studies were carried out in the quan- 
tum regime of multiple scattering using nonclassical light 
sources: spatial quantum correlations in the multiply 
scattered light were observed using squeezed light. The 
correlation function was determined by analyzing to- 
tal transmission quantum noise measurements result- 
ing into an average over all spatial and angular direc- 
tions [12]. Recently, Peeters et al. observed two-photon 
speckle patterns by sending two spatially entangled pho- 
tons through diffusers [16] and spatial correlations were 
observed in photon coincidence measurements exploiting 
surface scatterer for a single configuration of disorder. 

Here, we access for the first time photon correla- 
tions of multiple scattered light in three-dimensional ran- 
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dom media directly using angular-resolved photon coin- 
cidence measurements between separate multiple scat- 
tering channels. The angular photon correlation func- 
tion is induced by multiple scattering of light and consti- 
tutes the multichannel analogue of the Hanbury Brown 
and Twiss experiment generalized to a multiple scatter- 
ing setting [17]. So far, such an extension of the two- 
channel correlation function has been actively studied in 
the realm of electron transport in mesoscopic conduc- 
tors [18]. We use a pseudothermal light source to inves- 
tigate the angular dependence and record the angle and 
time correlations of the multiply scattered photons. Af- 
ter a configurational average over realizations of disorder, 
we find that the angular photon correlation function is 
dependent on the photon statistics of the light source and 
on classical intensity correlations induced by the multi- 
ple scattering medium. Our experimental results demon- 
strate the validity of the general quantum formalism of 
multiple scattering based on a continuous mode theory 
in the time domain. We underline that our measurement 
scheme is not only restricted to the regime where a clas- 
sical description of multiple scattering suffices but could 
pave the way to study quantum interference phenomena, 
as theoretically proposed in Ref. 9. Furthermore, the 
investigation of photon correlations in speckle patterns 
could be of importance for the field of ghost imaging [19- 
21]. 

The dimensionless angular photon correlation function 
between two angular directions bo and bi (sec Fig 1), is 
defined as 



(1) 



where the classical ensemble average of the stochastic 
process of multiple scattering is denoted by the bars. 
Performing the ensemble average enables extracting sta- 
tistical correlations that are hidden in the fluctuations 
of the speckle pattern. The normal ordered correlation 
function in the numerator can be measured by record- 
ing the coincidence counts between two detectors within 
the time interval At [22]. The dimensionless mean pho- 
ton number is obtained by integrating the photon flux 
operator over the measurement time At: n{t, At) ~ 
dt' a^{t')a{t') [23]. We furthermore consider a sta- 
tionary light source whose statistical fluctuations do not 
change in time [24]. The correlation function there- 
fore only depends on the time difference t = t' — t oi 
the measurement between bo and bi and the time ar- 
gument t can be neglected. In the measurements, a 
single photon counting detector records the number of 
photons within a time interval At. Both (??(Ai)) as 
well as the variance in the photon number fluctuations, 
An'^{At) = {n^{At)) - {n{At))^ , are obtained by count- 
ing photons in many time intervals of length At. The 
ratio, F{At) = An^ {At) / [n{At)) , is called the Fano 
factor and gauges the magnitude of the photon fluctu- 
ations. F{At) > 1 corresponds to light sources with 
super-Poissonian photon statistics that are of concern 
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FIG. 1: a. Sketch of the experimental setup. L: lens, G: 
rotating ground glass plate, I: iris, P: polarizer, F: 10 nm 
interference filter, AO: angle between detector Di and D2; a: 
incident light channel; 60, 61: exit light channels. Optionally 
a neutral density filter (OD) is inserted in front of detector 
Di attenuating the light beam. The relative angle between 
D2 and the optical axis of the incident light beam is 4.5°. To 
probe different realizations of disorder the sample is displaced, 
b. Recorded photon counting distribution, P{n{At)), of the 
light source on detector Di after the neutral density filter and 
without the sample using a time interval of At = 266 fis to 
bin the measured photons. The photon counts n{At) have 
been measured TV = 11.000 times, c. Same measurement as 
in b, but using a time interval of At = 100 /iS. 



here, while F{At) < 1 is found for purely nonclassical 
light sources exhibiting sub-Poissonian photon statistics. 

Multiple scattering of super-Poissonian light is pre- 
dicted to lead to photon correlations between different 
angular directions of the output beams. Having an inci- 
dent light wave in direction a, we derive for r = [25] 
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The first term on the right side is related to the photon 
statistics of the light source and is the equivalent to the 
second order coherence function, g^'^\T = 0) [24], while 
the latter factor. 
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is due to classical intensity correlations [2]. In the dif- 
fusive regime classical intensity correlations vanish and 

^bb' ~ 1- ^^'^ special case of a thermal light source 
Eq. (2) simplifies to [24] 



ClUir,At) = [l + exp(-^(r/r,)2)] x C^^^, (4) 

with Tc being the coherence time. In the experiment At 
can be chosen in such a way that the Fano factor of the 
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FIG. 2: Classical intensity correlation function, C^^^^ ('^Q) 
versus AO for the speckle pattern shown in the inset. The 
full width at half maximum represents the average width of 
the speckles, a. The speckle pattern is obtained by rotating 
detector Di behind the sample. = 0° corresponds to the 
configuration where Di is positioned on the optical axis of the 

incident light source. The uncertainty in C^^^^ (^Q) does not 
depend on the angle and yields ±0.4. 

light source, F{At), corresponds to thermal light, i.e. 
F{At) = {h{At)) + 1. In contrast photons of a coher- 
ent light source are uncorrelated {F{At) = I) and the 
angular photon correlation function equals the classical 
intensity correlation function, cf. Eq. (2). 

The experimental setup is displayed in Fig. I. As a 
light source we use a continuous wave TiiSapphire laser 
(A = 780 nm) that is focused onto a ground glass plate. 
Super-Poissonian photon statistics is obtained by super- 
imposing coherent beams with random amplitudes and 
phases. For that purpose the ground glass plate is slowly 
rotated and only a fraction of the transmitted light is 
collected using an iris with an aperture that is smaller 
than the average speckle size of the scattered light gener- 
ated by the ground glass plate [26] . The iris ensures that 
only the part of the light is collected for which wave vec- 
tors exhibit maximum momentum correlation [27]. The 
light beam is then coUimated and directed onto the front 
surface of a multiple scattering medium consisting of tita- 
nium dioxide nano-particles (thickness, L = 6.3 ±0.2 fim, 
transport mean free path, £ = 0.9 ±0.1 fim [28]). Two 
single photon counting detectors (Di, D2) are positioned 
behind the sample to record photon coincidence counts 
between multiply scattered photons at different angles. 
Detector Di can be rotated to vary the angle A8 be- 
tween both detectors. We avoid the contribution from 
ballistic propagation of light through the multiple scat- 
tering medium by only collecting light polarized perpen- 
dicular to the incident light polarization. 

We ensure to probe different speckles, i.e. different out- 
put modes 60 and bi , by determining the average speckle 
width, a. A typical intensity speckle pattern of multiply 
scattered light is shown in the inset of Fig. 2, recorded 
by rotating the detector Di. From the angular-resolved 



speckle pattern we obtain the classical intensity correla- 

(c) 

tion function, C^^5^(A0) [3, 29] by calculating the auto 
correlation function of the speckle pattern versus angle, 
AQ, plotted in Fig. 2. The full width at half maximum of 
the classical intensity correlation function corresponds to 
the average speckle width, which is a = 1.5°. The corre- 
lation function approaches C^^^"*^ (^0) = 1 for A8 « 1.25 
showing that different speckles are uncorrelated when 
considering only intensity measurements as opposed to 
the correlation measurements reported below. The in- 
crease in the correlation function for large angles origi- 
nates from the limited number of investigated speckles 
and could be suppressed by using better angular filtering 
with a smaller iris in front of the detector or by prob- 
ing a larger angular range. The apertures in front of the 
photon counting detectors Di and D2 are adjusted to be 
smaller than a while in the following experiments, the 
angle between both detectors, AQ, is chosen to be much 
larger than a. Thereby independent output modes are 
probed and the classical intensity correlations within a 
single speckle (see Fig. 2) do not contribute to the angu- 
lar photon correlation function. 

In order to characterize the properties of the light 
source, we measure the photon statistics, i.e., we cre- 
ate a histogram from the number of photons, n(At), 
that are detected within At. To this end, the sam- 
ple is removed and Di is positioned in the light beam 
emerging from the source after attenuation with a neu- 
tral density filter. Figure lb displays the measured 
supcr-Poissonian photon statistics for At ~ 266 fis yield- 
ing FoD{At) = (noD(Ai)) = 29. The correspond- 
ing Fano factor of the light source without attenua- 
tion is calculated using the separately measured opti- 
cal density 0D= 2 x 10~^ of the filter and we obtain 
F{At) = (FoD(Af) - 1)/0D ± 1 = 1.4 X 10^ equal to 
the mean number of photons, {n{At)) — (noD(At))/OD. 
Hence, the light source incident on the multiple scat- 
tering medium exhibits the properties of thermal light 
and is called a Gaussian radiation source [26]. By vary- 
ing At, the photon counting distribution changes. Fig- 
ure Ic shows for example the measured (attenuated) pho- 
ton statistics for At ~ lOO/is with a corresponding Fano 
factor of Foi){At) ~ 17 and mean number of photons 
of (noD(At)) = 10. By removing the ground glass plate 
in the incident light beam, we measure FoD(At) = 1, 
independent of At, reflecting the coherent state of the 
continuous wave Ti:Sapphire laser. 

The correlation between photons that exit the multi- 
ple scattering medium under a relative angular difference 
of AQ is obtained by recording simultaneously the mean 
number of photons within At in each direction, (fit, (At)). 
In addition we measure the photon coincidence counts, 
{h{,g{At)hi,-^{At)) , within the time-interval At between 
detector Di and D2, see Eq. (1). The ensemble average 
is achieved after repeating the measurement at 200 dif- 
ferent sample positions (indicated by the arrows at the 
sample in Fig. 1) corresponding to different realizations of 
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FIG. 3: Photon correlation function C^^^^ (0, At) (gray circles) 
depending on {h{At)) and F{At) (top axis) of the light source 
that is obtained by varying the time interval At while slowly 
rotating the ground glass plate, see Fig. la. For a coherent 

state C^jj^(0, At) is independent of (n(At)) (black squares). 
The angle between the detectors is chosen to be AS = 9° and 
the theoretical predictions contain no free fitting parameters 
(solid curves). 

disorder. The experimental configuration further allows 

(c) 

us to measure the classical intensity correlations, C^gti^j 
independently, as explained above. These measurements 
are fundamentally different from the photon coincidence 
measurements sine we record the individual ensemble- 
averaged mean number of photons (fib (At)) on each de- 
tector and the joint ensemble average (nto {^t)){f^bi i^t)) 
on both detectors. 

First, we study the angular photon correlation func- 
tion, C^^^(0, At), depending on the average number of 
photons and Fano factor of the incident light field by 
varying At, as shown in Fig. 3. For light sources with 
super-Poissonian probability distributions {F{At) > 1) 
we measure a decrease in the strength of the photon cor- 
relation function with increasing number of photons and 
Fano factor (gray data points). On the contrary, we ob- 
serve that for a coherent light source, C^j^(0, At) is in- 
dependent on the incident number of photons. In order 
to compare the experimental data with theory we mea- 
sure also the classical intensity correlations induced by 
the multiple scattering medium. Afterwards, the Fano 
factor of the light source is determined with detector Di 
as a function of At by removing the sample. The ex- 
perimental data are found to be in good agreement with 
Eq. (2) without any adjustable fitting parameters. The 
slightly smaller predicted values of the angular photon 
correlations (gray line in Fig. 3) are attributed to varia- 
tions in the rotation speed of the ground glass plate over 
time, which influence the Fano factor. For a very low 
number of incident photons we observe a saturation at 

C^^^(0, At) w 2.9 corresponding to very short time in- 
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FIG. 4: a. Angular dependence of C^^j,^ (0, At) for a coherent 
light source (black squares) and a Gaussian radiation source 
(gray circles). Detector Di is rotated behind the sample to 
change AO. The straight lines are the theoretical predictions 
assuming that the classical intensity correlations equal unity, 
b. Measured temporal dependence of the angular photon cor- 
relation function for a Gaussian radiation source (gray cir- 
cles) and a coherent light source (black squares). Each data 
point represents an individual measurement using A0 = 9° 
and At = 266 /xs. For visibility the uncertainties of the data 
points are not plotted which are ±0.03 for the coherent light 
source and ±0.07 for the Gaussian radiation source. The gray 
curve is a fit to the data with the coherence time as a free pa- 
rameter. 



tervals At. This effect we attribute to the properties of 
the light source determining its photon statistics. The 
phase of an otherwise coherent light source is scrambled 
by the rotating glass plate imposing an upper limit on the 
ratio (F(At) — l)/(n(At)) and the correlations, respec- 
tively, for times At smaller than the time scale on which 
the phase distortion occurs. These results demonstrate 
that the photon fluctuations of multiple scattered light 
depend sensitively on the investigated timescales, even 
though samples with static disorder are probed. Thus, 
the strength of the angular photon correlation function 
can be controlled by varying the bin time At in the mea- 
surements. 

Figure 4a plots the main result of this report, i.e. the 
angular photon correlations as a function of the angle 
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A8 between detector Di and D2 using a coherent and 
a Gaussian radiation source with At = 266 /iS. Here 
detector Di is rotated behind the sample, as shown in 
Fig. la. For each detector position we recorded for 200 
different reaUzations of disorder the individual as well 
as joint photon counting distributions (see Fig. lb). We 
neither observe for the coherent light source nor for the 
Gaussian radiation source a dependence of the angular 
photon correlation function on the angle. This holds in 
the diffusive regime of multiple scattering, where the clas- 
sical intensity correlation function does not vary with an- 

(C) 

gle, C^obi " ^- ^'^^ ^ source with super- Poissonian 
photon statistics the photon correlation function differs 
from unity demonstrating that different far-field speckles 
are strongly correlated and the ensemble-averaged cor- 
relation does not decrease as the angle increases. We 
underline that the relative angle, A8, between the de- 
tectors is much larger than the speckle width, a, i.e. we 
probe different speckles that are classically uncorrelated. 
These new correlations do not depend on the angle be- 
tween the speckles and are therefore infinite in range, in 
agreement with theoretical predictions (straight lines in 
Fig. 4a). Thereby, we show that a multiple scattering 
medium translates the temporal photon correlations of 
the light source into angular correlations connecting si- 
multaneously all input and output modes to each other 
- a property that has been suggested as a way to gener- 
ate quantum entanglement in a disordered medium [7, 9] . 

The error bars in C^^^(0, Ai) arise mainly from uncer- 
tainties in the classical intensity correlations caused by 
a finite number of disorder realizations for the ensem- 
ble average. In previous work on with squeezed light 
the photon correlation function was extracted from total 
transmission quantum noise measurements [12], and not 
through direct photon counting as employed here. Wc 
note that the technique of measuring the quantum noise 
of the total transmission to extract the photon correla- 
tion function only applies in a diffusive regime, while the 
direct detection scheme developed here can be applied 
also in mesoscopic and localized regimes where interest- 
ing new quantum correlations are predicted [5]. Fur- 
thermore, direct photon counting measurements would 
be required for testing recent predictions that quantum 
interference can survive multiple scattering even after en- 
semble averaging over configurations of disorder [9] . 

Our experiments on the photon correlations after mul- 
tiple scattering are potentially relevant to the field of 
ghost imaging, where the second order coherence be- 
tween two separate, but spatially correlated light beams, 
can be used to spatially resolve an object in one of the 
beams [20]. A light beam emerging from a pseudother- 
mal light source, similar to the one used here but without 
an iris behind the ground glass plate (see Fig. la), is di- 
vided in two optical paths using a beam splitter [27]. In 
one of the paths the object is placed in front of a large 
bucket detector and the second light beam is imaged with 
a lens onto a photo detector with a small aperture. Since 
this light source is chaotic, only photons within a single 



speckle spot of the light exiting the ground glass plate 
are spatially correlated, resulting in point-to-point cor- 
respondence of the photon correlations between object 
plane and image plane [20]. By scanning the detector 
across the image plane, the object can be resolved from 
the spatial dependence of the photon correlations. Our 
experiments are distinct to this scheme since the pho- 
ton correlation function observed here extends over the 
entire speckle pattern that is generated by multiple scat- 
tering, and the implemented iris ensures that only pho- 
tons of the light source with similar wave vectors are 
selected. The infinite range of the angular photon cor- 
relation thus originates from the photon statistics of the 
incident light beam distributed over the available degrees 
of freedom due to the multiple scattering [5, 18]. Thus, if 
multiple scattering occurs after the pseudothcrmal light 
source, the point-to-point mapping of the second order 
coherence is lost and cannot be exploited for ghost imag- 
ing. A second approach to ghost imaging does not utilize 
photon fiuctuations of the light source but uses a multi- 
ple scattering medium in front of the beam splitter. In 
this case the generated intensity speckle pattern is col- 
limatcd and split in two light beams, establishing now 
a point-to-point correspondence of the classical light in- 
tensity between object and image plane. After perform- 
ing an ensemble average over many different speckle pat- 
terns, the resultant classical spatial intensity correlation 
function resolves the object [21]. These two fundamen- 
tally different ghost imaging approaches utilize the pho- 
ton correlations of the light source and intensity correla- 
tions induced by multiple scattering, respectively. Here, 
we use a pseudothcrmal light source to probe angular 
photon correlations that are induced by a multiples scat- 
tering medium and study its dependence on the photon 
fluctuations of the incident quantum state of light, cf. 
Fig. 3. While ghost imaging schemes so far exploit cor- 
relations between two light beams, multiple scattering 
would in general allow to induce correlations between 
output modes obtained by mixing N independent quan- 
tum states. 

In order to investigate the temporal dependence of the 
angular photon correlations induced by a multiple scat- 
tering medium, a time delay, t, is induced between de- 
tector Di and D2. Figure 4b presents the angular pho- 
ton correlation function depending on r using a Gaussian 
radiation source (gray data points). As the time differ- 
ence increases, a clear decay of (t. At) is observed. 
Fitting the experimental results with theory (Eq. (4)), 
we determine the coherence time to be ~ 750 fis. 
The coherence time describes the dephasing of a light 
source, i.e. the average time interval between phase dis- 
tortion fluctuations. Only coherence times of the light 
source can be resolved that are larger than At, i.e., the 
experimentally measured photon statistic does not con- 
tain information about photon fluctuations on shorter 
timescales. As Tc > At we resolve the coherence time 
of our Gaussian radiation source. The time response 
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of C^gti (''"' originates from the temporal bunching 
of photons in the thermal state. For comparison, the 
temporal dependence of the angular photon correlation 
function is investigated for a coherent light source, too 
(black data points, Fig. 4b), where as expected the angu- 
lar photon correlation function is independent of r and 
equals unity. 

In conclusion, we demonstrated experimentally the 
first direct measurement of angular photon correlations 
that are induced by multiple scattering of light. The time 
correlations of the multiply scattered photons at different 
angles were found to be dependent on the photon statis- 
tics of the incident light source. In excellent agreement 
with the quantum theory of multiple scattering, we ob- 



served that the photon correlations are infinite in range 
in the diffusive regime. Combining the quantum aspects 
of multiple scattering with phase shaping techniques [30] 
could open a new route to vary the quantum state of 
the multiple scattered light. Our angular-resolved exper- 
imental technique is crucial to study experimentally pre- 
dicted quantum interference and quantum entanglement 
in three-dimensional multiple scattering media [7-9]. 
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